Momentum-space analysis of multipartite entanglement at quantum phase transitions 
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We investigate entanglement properties at quantum phase transitions of an integrable extended Hubbard 
model in the momentum space representation. Two elementary subsystems are recognized: the single mode 
of an electron, and the pair of modes (electrons coupled through the ^-pairing mechanism). We first detect 
the two/multi-partite nature of each quantum phase transition by a comparative study of the singularities of Von 
Neumann entropy and quantum mutual information. We establish the existing relations between the correlations 
in the momentum representation and those exhibited in the complementary picture: the direct lattice representa- 
tion. The presence of multipartite entanglement is then investigated in detail through the Q-measure, namely a 
generalization of the Meyer- Wallach measure of entanglement. Such a measure becomes increasingly sensitive 
to correlations of a multipartite nature increasing the size of the reduced density matrix. In momentum space, 
we succeed in obtaining the latter for our system at arbitrary size and we relate its behaviour to the nature of the 
various QPTs. 

PACS numbers: 03.65.Ud, 71.10.Fd 



I. INTRODUCTION 

The idea of studying and characterizing Quantum Phase 
Transitions (QPTs) by means of tools and concepts borrowed 
from the quantum information theory has provided in the past 
years a significant amount of interesting results'. In partic- 
ular, the notion of entanglement and its measures is of cen- 
tral importance within this framework. Indeed, for several 
and different models in condensed matter physics it has been 
shown that partial derivatives of appropriate measures of en- 
tanglement with respect to the field driving the QPT keep 
track of the underlying critical phenomenon by developing a 
singularity 2,3 ! 4 ' 5 ' 6 ' 7,9,10 - 11 ' 12 ! 13 ! 14 . Moreover, from an entangle- 
ment study of QPTs one can gain further information such as 
the order or the scaling exponents as well as a deeper insight 
on the structure of the correlations involved. For instance in 10 
it has been argued how, for an extended Hubbard model, a 
comparison of the Von Neumann entropy -a measure of quan- 
tum correlation- of a single site with the quantum mutual in- 
formation, measuring the total correlations between two sites, 
allows one to infer which kind of correlation (two-point or 
multipartite) plays the prominent role at a given quantum tran- 
sition. Such a scheme has been tested i n 14 ! 15 . In particular, 
ini 4 - the two-point QPTs have been investigated by means of 
punctual measures of two-point quantum correlations. As for 
the multipartite ones, the analysis becomes harder. In fact, the 
evaluation of the known direct measures of multipartite cor- 
relations requires either some optimization process that turns 
out to be in practice unfeasible, or the knowledge of the re- 
duced density matrix relative to arbitrarily large subsystems. 

Entanglement is in general a notion that depends on the 
choice of the subsystems whose correlations one is interested 
in analyzing. In particular, the same Hilbert space Ti can 
support different tensor product structures and thus different 
choices of subsystems are possible 23 . In describing many- 
body condensed matter systems one is naturally led to con- 
sider at least two different structures: the one induced by the 
real-space topology of the system under consideration, i.e., the 



direct lattice picture, and the one corresponding to the mo- 
mentum space, i.e., the reciprocal lattice picture. Although 
the literature on entanglement in condensed matter has privi- 
leged the first description, one expects that the information ob- 
tained by the analysis of entanglement properties of elemen- 
tary subsystems of the direct lattice (subsets of single sites of 
the ambient lattice) can be fruitfully integrated with that ob- 
tained by considering as elementary subsystems the momenta 
modes of the reciprocal lattice. 

In this paper we provide the momentum space description 
of the ground state entanglement properties of the extended 
Hubbard model discussed i n 1( ¥ 4 ' 16 . The study in momentum 
space allows for a different and rich picture of the structure of 
correlations and of their relations with QPTs. Indeed, the rel- 
evance itself of multipartite (multi-mode) and two-point (two- 
mode) correlations at a given QPT can significantly vary when 
the above change of perspective is performed. Furthermore, 
the comparison of the (quantum) correlations in the momen- 
tum and direct lattice picture turns out to be a fruitful method 
of analysis. Indeed, this procedure allows one to enlighten 
the origin of the off-diagonal long-range order (ODLRO) dis- 
played by the ground state for appropriate values of the param- 
eters: this is directly related to the value of negativity between 
modes k and -bi2. 

Moreover, the possibility of determining the reduced den- 
sity matrices of arbitrary subsystems composed by sets of mo- 
mentum modes, allows one to overcome the above mentioned 
difficulties in measuring directly multipartite correlations. 
We exploit such possibility to compute the generalization of 
the Meyer- Wallach measure of multipartite entanglemen t 20 ' 22 . 
We thus show that by increasing the subsystem size the partial 
derivative of Q develops a singularity in correspondence with 
precisely those QPTs at which the dominant correlations were 
conjectured to be of multipartite nature. 

The paper is organized as follows. In section|II]we describe 
the phase diagram of the model under investigation and we 
introduce the momentum space description its ground state in 
the various phases. In section HIT1 we first discuss the general 
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structure of the correlations in the ground state in the whole 
phase diagram; we related the picture emerging in the mo- 
mentum space with the one previously obtained in the direct 
lattice picture; we discuss the role of the correlations at QPTs 
and we evaluate and discuss the Meyer- Wallach measure of 
multipartite entanglement. In section [TV] we summarize our 
conclusions. 



reads (n s — N s /L,n d — N d /L): 

e(n s ,n<i) = — sin (wn s ) + una ■ (3) 

The minimization equations can be solved, thus obtaining 
the ground-state phase diagram shown in Fig. Q]in the n-u 
plane (with n = N/L average per-site occupation number). 
Three distinct regimes are recognized, depending on n s , rid- 



II. THE BOND-CHARGE EXTENDED HUBBARD MODEL 
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The model we investigate is the bond-charge extended Hub- 
bard model for a particular value of the bond-charge interac- 
tion at which the model becomes integrable 16 . In units of the 
nearest-neighbors hopping amplitude, its Hamiltonian reads: 
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H = - ^2 [1 - ( n is + n js )]c\ a c oa + u^n^n ll , (1) 
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where cj CT and c ia are fermionic creation and annihilation op- 
erators on a one-dimensional periodic chain of length L; a = 
{|, 1} is the spin label, a denotes its opposite, n = c] c - 
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is the spin-cr electron charge, and (i, j) stands for neighbor- 
ing sites on the chain; u is dimensionless on-site Coulomb 
repulsion. 

The two following feature make Hamiltonian (Q~|i 
integrable 16 : i) the number of doubly occupied sites is 
a conserved quantity; ii) the role of spin orientation is 
irrelevant. For an open chain, any sequence of spins in the 
chain cannot be altered by Hbc, while for periodic boundary 
conditions, only the sequences of spins related by a cyclic 
permutation can be obtained. In particular, the ground state 
turns out to be degenerate with the fully polarized state. 



A. Eigenstates, eigenvalues and ground-state phase diagram 

The physics of the system described by H is basically that 
of N s — (N s ) spinless fermions moving in a background of 
L — N s bosons, of which Nd = (N d ) are pairs of particles 
and the remaining are empty sites. Both N s and Nd are con- 
served quantities (and determine the total number of electrons 
N = N s + 2Nd) and the following ansatz for the form of the 
eigenstates can be formulated in momentum space: 
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\ip(N s ,N d )) 



■ Af(tf) N *a 



t 



(2) 



Here TV = 



(L — N s — N d )\/{L — N s )\N d \ 



1/2 



is the nor- 



malization factor and a\. a = J^i ' 

is the so-called eta operator, its action 
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being to create a 0-momentum couple of fermions delocalized 
over the reciprocal lattice. 

The actual ground state \ipcs(N s , Nd)) is chosen among 
the eigenstates (O by requiring that N s and N d minimize the 
corresponding eigenvalue, which in the thermodynamic limit 
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FIG. 1: fc-space representation of the ground-state phase diagram 
of H in the n-u plane. Dotted lines represent some cases of iso- 
correlation curves, see sec. IIIIBI 

For sufficiently high u, n d = characterizes the conducting 
region /, which, for n = 1 becomes insulating (phase 
IV). In this case n s — n. Keeping n fixed and decreasing 
u < u c (n) — — 4cos(7rn) we pass to region II, where 
the pairing mechanism becomes convenient and the energy 
minimum is realized for rid ^ 0. By further decreasing u 
we enter in region III, where pairs are energetically favored 
with respect to single fermions; the energy minimum is thus 
obtained imposing the condition n s = (i.e. rid — n/2). Al- 
though both phases I and II are characterized by a conducting 
behavior, they are qualitatively different since only in phase 
II the occurrence of off-diagonal long-range order (ODLRO) 
- which also survive in phase III - shows up 14 . 

The standard approach to the study of (quantum) phase 
transitions consists in investigating the behavior of the 
(ground state) energy density at critical points. Within this 
framework, all QPTs here studied are second order ones: they 
are signaled by non analyticities of the second order partial 
derivatives of the ground-state density energy. When u = ±4 
two transitions take place: II^IV and II— > III. Here the signif- 
icant derivative is the one with respect to the field u. Indeed, 
it is the second derivative that diverges at the critical points 
u r . = ±4: 



9lE gs (u,0) 
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(4) 



At the transition II— >I the QPTs are again signaled by the 
non analyticities of the second order partial derivatives, but 
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-at variance with the previous case- these simply exhibit 
a discontinuity and do not diverge at the transition points 

u c = —4 cos(7m): 

lim d 2 u E gs (u,n) = ? - ^ = = 
u^uj OTTy/1 — cos' 1 (mr) 

= lim _dlE gs (u,n) , 

U — *U C 

lim d^E gs (u,n) = 7^ 2tt cos (n7r) = 

n — >n~£ 

= lim_ d^E gs (u,n) . 

u — m c 

(5) 

As already mentioned, the QPTs above described were 
studied inii in terms of (quantum) correlations between dif- 
ferent subsystems in the direct lattice picture. There the natu- 
ral subsystems were the single/pairs of sites of the lattice. In 
this picture, the principal measures of correlations taken into 
account were: i) the Von Neumann entropy of the generic sin- 
gle site I, Si, that measures its total quantum correlations with 
the rest of the chain; it) the mutual information between two 
sites I, m, Zj >m , that measures the total (quantum and classi- 
cal) correlations between the two sites; Hi) the negativity 17 
(and concurrence 18 ) between the two sites l,m, A/j, m , that 
measures only the quantum contribute to the correlations be- 
tween the two sites. Remarkably, the study of the derivatives 
of these measures of correlations allowed to distinguish the 
various QPTs in terms of relative relevance of the two-point 
(i.e., two-site) and multipartite correlations. In particular, the 
divergence (HJl exhibited by the energy density is the same 
shown by the first derivative of both quantum mutual infor- 
mation and Von Neumann entropy at both transitions. This 
is to understand as a signature of the relevance of two-point 
entanglement at the transitions, relevance that was confirmed 
by the divergence of the entanglement ratio that, indeed, orig- 
inates from the spreading of the two-point quantum correla- 
tions along the whole chain. 

In correspondence with transition II— »I -signaled by the 
discontinuities of the energy density (0- while the partial 
derivative of the Von Neumann entropy still exhibits a loga- 
rithmic divergence, the partial derivative of the quantum mu- 
tual information is simply discontinuous. This is to under- 
stand as a manifestation of the relevance of multipartite en- 
tanglement at such transition. 



B. fc-space representation 

The main aim of this paper is to extend the study of how the 
QPTs introduced in the previous sections can be described in 
terms of correlations in the fc-space picture. The fc-space rep- 
resentation of the eta-paring states with Nn fixed and N s = 
have been thoroughly analyzed in 1 -. Here we first give a brief 
description of the main ingredients needed for the analysis; 
second we see how the picture can be extended to region II in 
which N s 0, i.e., the ground state is filled at the same time 
by single fermions and pairs of particles. 



As previously mentioned, each single eta operator in the the 
fc-space representation can be written as: 

3 3 

where fj k . = 0—k-t a k-l m ^ a ! cr' afc j cr are f erm i° n i c opera- 
tors that create/annihilate a particle with given spin a kj in one 
of the L momentum modes kj = 2wj /L whose local basis 
is B kj = {\0)k f ,\ T)fy,| U)fej}- First, let us notice 

that the different regions in Fig. Q] correspond through (O to 
different realizations of the ground states in fc-space. In partic- 
ular, in region /, at variance with the direct lattice picture, the 
ground state, being constituted by N s electron with definite 
momentum, is completely factorized: 

KW) = <=iK> ■ (7) 

In region II the ground state (O is built by occupying with 
N s fermionic particles (with fixed spin a kj ) the modes of the 
reciprocal lattice 

kj € K k < kNs = {K min = -ir(N s - 1)/L, . . . , 

^max - tt(N s - 1)/L} . (8) 

Since in region II part of the momentum space is occupied by 
N s single particles, the derealization implemented by the eta 
operators can only take place over a subset of momenta: 

K k >k Ns = {ttN s /L,...,tt} 

U{-tt + 2-k/L, -2ir{N s /2 + l)/L}. 

The number of distinct values of momenta available for the 
creation of pairs of particles, i.e., the order of K k >k N , is 
L — N s . This fact can be represented by observing that the 
partition of set of possible fc-modes in K k < kN , ^K k >k Ns cor- 
responds to a bipartition of the fc-state space: Ti = Tik<k Ns ® 
T~Lk>k N , ■ Considering the vacuum 

|vac) w = |vac) Hfc < fcjVs ® |vac) Hfc>fcjVs (9) 

the ground state in (O can be written as 

\1>n(N„N d )) = Af(v f ) Ni al <T ...a{ N Jx aC ) H 

= (rjt)^|vac> Wjb> ^J^(JV.)>. 

^From the point of view of the derealization implemented 
by the eta-operators, region III represents the limiting case 
where the 2V<j pairs can be created over the full momentum 
space: 

|feW)iAr(?,t)^|vac)w (10) 

The above discussion about region II implies that the only 
modes that are non-trivially correlated are those belonging to 
T~(-k>k N ,, while the modes in Tt k <k N ^ are simply in a tensor 
product state where each mode is singly occupied by a spin-er 
particle. The characterization of region II and of the transi- 
tions II — can thus be done in terms of the correlations 
between the modes over which the pairs of particles created 
by the eta-operators are delocalized. The influence of the sin- 
gle particles in this picture is to reduce of the space available 
for the derealization and this reduction has an important role 
in determining many properties of the state correlations. 
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HI. RESULTS IN MOMENTUM REPRESENTATION 

We now generalize the results obtained inl 9 . for the eta- 
pairing states with N s = to the generic states (f2]) where 
the eta operators act on a reference state characterized by the 
presence of single particles with a given momentum. The ex- 
tension allows one to investigate, in fc-space, the behavior of 
the correlations at the QPTs exhibited by the model ([I]). 



A. Local structure of correlations in k-space 

A first significant result in the analysis of the structure of 
correlations is the acknowledgement that both the generic sin- 
gle mode kj and the paired modes (—kj, kj) are fundamental 
building blocks. The basic objects that one has to evaluate 
in order to derive the appropriate measures of correlations are 
the reduced density matrices of the given subsystems: i) p k . 
for the single mode; ii) p ki , kj f° r two generic modes; Hi) 

P(-ki,ki),(-k -,k ) f° r two P a ii" s °f modes. 
The general considerations that allow one to determine the 
structure of these density matrices are the following. Since 
the state (rf) a \\ac)-n k>kN is built by the repeated applica- 
tion of operators rj\ , each time the — kj mode is occupied, 
the kj one must be occupied correspondingly. Therefore, for 
example, the only non-vanishing elements of p k . are the di- 
agonal ones. As far as p k% . kj (p(_ kitki) ^_ k .^) is concerned, 
the only d-modes (d = 2, d = 4 respectively) projectors P^ d 
that have non-vanishing expectation values are those preserv- 
ing the number u d = 0, . . . , d of pairs of particles (a k . , . 
For a given v& the expectation values are all equal and in the 
TDL they read 



(i**) = a Vd (l - a) d - Ud , where a = 



n d 



1 - n. 



(ID 



It turns out that these calculations can be generalized to an ar- 
bitrary number of both single and paired modes, one can thus 
explicitly find the diagonal form of the reduced density matrix 
of any subsystem. In the TDL, die latter and consequently all 
measures of correlations depend exclusively on the parameter 
a. 

Some basic considerations are here in order. As previously 
observed, the presence of the N s single particles with fixed 
spin a k . (kj < fcjv s ) in region II has the effect of reducing the 
available modes over which the eta-operators can delocalize 
the pairs of particles they create. This effect can be substan- 
tiated by considering the average number of particles with a 
given momentum k 3 ■: (n k .) = (nkA + n k . 1) . The latter is 
simply constant and equal to one for all kj € K k < kN ■ For 
K k >k Ns we have that (n kj ) = 2 (n kja ) = 2N d J(L-N s ) = 
2n d /(l — n s ) = (n — n s )/(l — n s ) = 2a. This means that, 
in general, the effect on correlations of the derealization im- 
plemented by the eta-operators depends only on the density of 
particles, i.e., it is related to the filling N — N s + 2N d and 
the number of available modes \K k>kN | = L — N s . 

We can now turn to the detailed description of the structure 
of correlations at a local level. In particular, the generic mode 



kj is quantum correlated with the rest of the system and, in 
the TDL, the Von Neumann entropy reads S kj = — 2 [a log a + 
(1— a) log (1 — a)]; furthermore kj has two-mode correlations 
only with the mode —kj'. 



S ki + 2a(l — a) ki = — kj 
k<i y— k j 

a(l — a)/3 ki = —kj 
h 7^ 



(12) 
(13) 



where I ki . kj (Af kitkj ) is the mutual information (negativity) 
between two generic modes (ki, kj). 

As for the paired modes (—kj, kj), each of them exhibits 
quantum correlations only with the rest of the system as a 
whole, indeed, while their Von Neumann entropy reads 



-2 (a log a + (1 - a) log(l - a)) + a(l - a) 



(14) 

the two-pairs correlations Ti_ k . k )j- ki , ki ) = 2a(l — a) [2 + 
a(l — a)(31og 2 3 — 5)] are just of classical nature since 



= 0. 



The above picture naturally suggests to investigate at QPTs 
the behavior of the correlations concerning the single mode 
and the paired modes. As already anticipated the expression 
of (n k ) is at the basis of the definition and the behaviour of 
( [TTirBi and in general of the correlations between any pair of 
subsystems (in the limit of large L). Therefore, the derivative 
of the measures of correlations at the various phase transitions 
can be expressed in terms of the derivatives of (n kj ), and this 
turns out to have important consequences for the discrimina- 
tion of the various QPTs. 

The above description of the network of correlations in fc- 
space can be fruitfully completed by analyzing the relations 
between such kind of correlations and some meaningful quan- 
tities that are defined in the direct lattice picture. In particu- 
lar one can see how ODLRO, which is defined in the direct 
lattice, is in fact directly related to negativity between paired 
modes, generalizing the relation established in' 9 . Indeed, in 
this model the ODLRO is defined as 



lim 

-j|-< 



(jiiVj} = n <i(l - n d) OC I, 



III 

hi 



(15) 



where 77J = c| jet . is the operator that creates a pair of parti- 
cles at the lattice site i; the the proportionality to the quantum 
mutual information between two sites i,j evaluated in region 
III i.e., X] a 1 , was found in 14 . Here we find that 



lim (ifiVj) = n d(l - n s - rid) oc J\f ki - ki (l - n s ) 2 . 

\i—j\—KX> \ I 

(16) 



which means that, even in region II i.e., where n s ^ 0, the 
negativity and the ODRLO are directly related and the latter 
relation reduces to the one found in 19 when evaluated in re- 
gion III (n 3 = 0). 

The relation between quantities defined in the different pic- 
tures can be further extended by considering the full expres- 
sion of Tfj i.e., the quantum mutual information between two 
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sites i,j in region II. We recall that in^i we derived the fol- 
lowing expression 



jii 



T 1 



rll 



(17) 



where 3-u_ji is the mutual information between the site i,j 
evaluated in region I (rid — 0, n — n s ) and X r 



ii 



lim 



*ao i s tne mum al information between two in- 



finitely distant sites in region II. This equation suggests that in 
region II one can distinguish the part of the direct lattice two- 
sites correlations i.e., Xh , which are determined by the frac- 
tion of n s single particles which are not correlated in fc-space, 
by another contribution, X~* t which takes into account, in di- 
rect lattice, the effect of the eta paring mechanism with which 
the rest of the fraction of particles are created . Here we can 
now fully understand this peculiar behaviour in terms of quan- 
tum correlations between paired modes. Indeed on one hand, 
by using the elements of the two-sites reduced density matrix 
which were analytically evaluated i n 10 ' 14 , one can show that 
the proportionality established in ( fT5b can be in some sense 
extended in region II: we have that in this region the ODLRO 
coincides with X^. On the other hand, thanks to ( fT6l ), this 
quantity can be directly linked to the A/fej.-fcj and this roots 
the origin of this part of the two sites correlations in the exist- 
ing correlations between paired momentum modes. 
It turns out that this is important relation between the (quan- 
tum) correlations pertaining the two complementary pictures, 
will also be relevant in the next section when we discuss the 
behaviour of the various correlation measures at the QPTs. 



B. Entanglement and quantum phase transitions 

In previous paper s 10 ' 14 we showed how, in the direct-lattice 
picture, the transition II— >I (see Fig. [U can be distinguished 
from transitions II — >III and and I— >IV on the basis of the 
different kind of correlations (two-site versus multipartite) 
among different subsystems. Such a distinction between these 
QPTs can be achieved in the fc-space picture only in a differ- 
ent way. In fact, we first note that at variance with the direct- 
lattice picture, here the insulating state which represents the 
ground-state at half-filling for u > 4 (insulating phase IV), as 
far as correlations are concerned, is not qualitatively distinct 
from the conducting state (O which occurs at arbitrary filling 
when u > — 4cos(-7m). Indeed, in fc-space the conducting 
state (0 is fully factorized for any filling n, so that entan- 
glement is identically zero at any n and the transition I— >IV 
cannot be detected from entanglement analysis. In the ambi- 
ent lattice the factorization property holds at half-filling only. 

In Table [I] we report the behavior at QPTs of the relevant 
partial derivatives of the studied quantities: Von Neumann en- 
tropy, quantum mutual information and negativity. 

We start by taking as elementary subsystem the generic sin- 
gle mode kj. At all the considered transitions we can rec- 
ognize that the two-mode correlations play a prominent role. 
Indeed, the correspondence of the divergencies of d x Sk and 




FIG. 2: (Color online) Plot of 9 U 5„ 



(straight blue line) and 



d u Z< 



■u.-L(-k i ,k i ),(-kj,k : j) 



(dashed red line) in region II. 




FIG. 3: (Color online) Plot of dnS-k^kj (straight blue line) and 
9„J(_ fciifci)i( _ fcj , fc .) (dashed red line) in region II. 



dxXkj-kj allows one to acknowledge the two-mode total cor- 
relations between kj and — kj as relevant for the transitions 
II— >I and II— TIL The distinction between the latter resides 
in the fact that only at transition II — >III d u Nk 3 -k 3 is singu- 
lar, suggesting that only in this case the two-mode quantum 
correlations between the modes — kj and kj play a prominent 
role. We note, however, that the two-mode character of the 
transition II— Til is of a different nature with respect to the 
one seen in the direct lattice picture, where we observed a di- 
verging entanglement ranged. 

This distinction in enforced by the relation between the neg- 
ativity and the ODLRO and the fact that the derivatives of 
these two quantities display the same divergence at the tran- 
sition II— >III. The comparison between these quantities and 
the mutual information between two arbitrary sites of the di- 
rect lattice discussed at the end of the previous section, allows 
once again to establish a link between the two-mode picture 
in fc-space and the two-site picture in direct lattice. Indeed, 
due to the relation (JTTJ the divergence of d u I^_., at the tran- 
sition II— > III, that allowed in 14 to describe this transition as a 
two-point one, is in part determined by the divergence of the 
d u Nk 3 ,-k 3 2- 

We now consider the elementary subsystem formed by 



the paired modes (—kj,kj). The vanishing negativity 
jV(-ki,ki),(-k-,k ) — shows that the two-pair correlations 
measured by i(-k ( ,ki),(-k-,k-) are always classical; this sug- 
gests the presence of just a multi-pair structure of quantum 
correlations. The divergencies of d x=u , n S-k .k show that 
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these correlations play an important, though different, role in 
both kind of transitions. The different role can be described at 
the level of two-pair correlations allowing for a classification 
of the QPTs into two classes: two-pair (transition II — ^III) and 
multi-pair (transition II— Indeed, as shown in table U and 
in figure [3] the singular behavior of d x Z(-kiki),(-k-k-) an d 
dxSkj-kj coincides only at transition II— > III. As in the direct 
lattice picture, we can then conclude that also in fc-space that 
the transition II — >III is characterized by the fact that correla- 
tions of two-point (two-pair) nature are the more sensible to 
the undergoing phase transition. On the contrary, the absence 
of a divergence of d x I^-^ki),{—k-k-) suggests that correla- 
tions of multipartite (multi-pair) nature are more relevant at 
the transition II— >I. In the next section we will provide a fur- 
ther confirmation of such statement by directly evaluating a 
measure of multipartite correlations. 



We now pass to the comparison of the singularities of the 
partial derivatives of the measures, as reported in TableQ] with 
those shown by the second partial derivatives of ground-state 
energy, (0, ©. Interestingly, at transition II >III all the com- 
puted quantities develop the same singularity of d^E gs (u, n). 
On the other hand, at transition II— >I d^ =u n E gs (u, n) ex- 
hibits a finite discontinuity and the same behavior is shown 
by both negativity and quantum mutual information, while the 
Von Neumann entropy is logarithmically divergent. 

To summarize, it seems that when the transitions are sig- 
naled by divergencies the main role is played by "two-point" 
correlations. Instead, in the case of simple finite discontinu- 
ities the dominant role is played by the multipartite quantum 
correlations. In both cases, the appropriate derivatives of Von 
Neumann entropy diverges. 





dxSkj , dxSkj.-kj ,dxlkj,-ki 


OxAfkj.-kj ,d x 1-kjkj,-kiki 




// -> I(x = It) 


log(u c - u) 


FD 


Multi 


II -* I(x = n) 


log n — n c \ 


FD 


Multi 


II -> III(x = u) 


— — u c 


— l/yju — U c 


Two 



TABLE I: Behavior of the evaluated partial derivatives at the various 
QPTs (left column). 'FD' stands for Finite Discontinuity, 'Multi' for 
Multi-pair, while 'Two' stands for Two-pair. 



We conclude this section by considering the half filling 
case. The latter has to be treated separately since no diver- 
gencies in the derivative of any measure of correlation can be 
observed at the transitions. This feature is a particular case 
of a more general phenomenon. As previously highlighted, 
in the TDL all measures of correlation depend on the average 
occupation number of each of the modes belonging to the set 
Kk>k N ■ Therefore the curves in the n, u plane defined by 
fixing (rikj) to a given constant value can be considered as 
iso-correlation curves (an example of them is given in figure 
[TJ. Along such curves the correlations of any accessible sub- 
system do not change. Consequently by crossing the boundary 
separating region II and III along these curves no divergence 
in the derivative of the correlations measures would be ob- 
served. 

In particular, the condition n = 1, i.e., the half filling case, 
defines a curve which is simply parallel to the u axis. Thus, 
while N s - that determines the number of available modes - 
and Nd can vary with u, \Tikj) is constant along this curve 
and all measures of correlations in region II are constant and 
fixed to the maximal achievable value (e.g., 5^. = 2). In or- 
der words, in the fc-space picture the transition between region 
II and HI, while it is well described by all measures of corre- 
lations for n ^ 1, at half filling (at variance with respect to 
the direct lattice case) and along the iso-correlation curves, it 
simply does not have any impact on the structure of the corre- 
lations of the state. 



C. Multipartite entanglement 

In order to evaluate directly the multipartite part of the cor- 
relations and their behavior at QPTs, we consider the gen- 
eralization of the Meyer- Wallach measure of multipartite en- 
tanglement Q(\tp)} 20 to multi-qudit states introduced by Scott 

• 22 

d m I 1 \ 

^ = ^i^-^E^J ■ as) 

Each multi index i = {ii, . . . , i m } identifies m elementary 
subsystems e.g., momentum modes in our case, that compose 
a given principal subsystem S m e.g., block of m modes. The 
latter is characterized by its reduced density matrix pq. The 
factor d m I (d m — 1) normalizes the measure to one; d is the 
dimension of the Hilbert space of each elementary subsys- 
tem and in the case of momentum modes it is equal to 4. 
Q m .d can be expressed as the average linear entropy of S m 
and thus quantifies the average entanglement between blocks 
of momentum modes and the rest of the system. The more 
the size of the blocks is increased, the more Q m _d can cap- 
ture and measure the multipartite correlations present in the 
system 22 '-! In order to obtain a proper measure of the latter 
it is thus necessary to be able to compute the reduced den- 
sity matrix of blocks of modes of arbitrary dimension. While 
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in the direct lattice this is possible only in region in k- 
space one can generalize the results obtained in 19 and extend 
them to region II. The extension is obtained again by recogniz- 
ing that, in /c-space, the role of the single particle with given 
spin <7fc. and momentum belonging to the set Kk j <k Ns is to 
reduce the momentum spectrum available for the derealiza- 
tion of the eta pairs to a subset of L — N s modes. One can 
thus derive the subsystems' density matrix spectrum and, by 
implementing the definition of ( TT~8T > one has for m = D: 



Q 



DA — 



1 - 



L-N t 
D 



-1 D 

I 

£> 2 =0 



£ /(£>2)Tr(p^)Tr(^ 2 ) 



A D - 1 



(19) 



where 



Tr^J 



2Di 

E 

M=0 
D 2 



i /2D 1 



AI 



(L-N S -2D^ -I 2 
[ Na — M ) 



V N d ) 



Tr(P73 2 ) = E 



,L-N S -D 2 
u 1 \ \ N d -a 



(L-N 3 
\ N d 



(20) 
(21) 



and 



nf = 2 /2 " 1 (^-^-2c 1 -2 J ) 



f(D 2 ) = 



Do) 



(22) 



is the number of equivalent partitions of D momentum modes 
into D\ = D — D2 single modes kj and D2 /2 paired of modes 
(— kj, kj). 

In figure [4] we plot the Q-measure (16) at quarter filling 
( L = 1000, N = 500) as a function of the Coulomb 
interaction u and for increasing values of the subsystem 
size D. The data confirm that despite the reduced portion 
of fc-space available for the construction of the eta pairs 
the multipartite entanglement induced by their presence is 
of fundamental importance in the whole region II. Indeed, 
for all values of u, we observe an increasing saturation of 
Qi,D to its maximal value when increasing D (for D > 4). 
Furthermore, for any dimension D of the block of modes, 
Qi : D is a monotonic decreasing function of u and thus, since 
nd = l/2(n — n s ) — l/2(n — arccos (— u/4)/tt), it is a 
monotonic increasing function of the number Nd of eta pairs. 



all correlations vanish, the multipartite ones change in a much 
faster way. Hence, the direct measure of multipartite correla- 
tions fully confirms the predictions of their role at QPTs, as 
obtained in sub-section IIIIBI We emphasize that a study at 
low values of D - like the one that would be feasible in the 
direct lattice picture - would have suggested the opposite con- 
clusion regarding the relevance of multipartite correlations at 
the transitions. 

We finally note that the determination of the reduced den- 
sity matrices po for blocks of modes of arbitrary dimension 
D allows for the evaluation also of the block entropy i.e., the 
Von Neumann entropy of pr>. This can be done by extend- 
ing the results obtained in 19 . This extension shows that even 
in the case analyzed in this paper the asymptotic behaviour 
of the block entropy depends on the type of modes included 
in the block. In particular, the relevant result is that when 
the block is composed by D/2 paired modes only, the block 
entropy grows logarithmically with D. This typical growth 
shows on one hand the presence of multipartite entanglement 
and on the other hand the peculiar multi-pair nature of this 
kind of entanglement. 




0,0 

-4,0 -3,5 -3,0 -2,5 -2,0 -1,5 -1,0 -0,5 0,0 

u 



FIG. 4: Plot of the Q-measure for a chain of length L = 1000 at 
quarter filling; D = 1, 2, 4, 8, 16, 32. 



IV. CONCLUSIONS 



The multipartite entanglement contribution at the various 
QPTs can be inferred by analyzing the behaviour of Q4.D at 
the critical points. Again, it is only for D > 4 that the correct 
qualitative behavior of Q4.D at large D is captured: while at 
the the transition II — the curves smoothly reach the (con- 
stant) value exhibited by the measure in region III, the plot 
clearly shows that d u Q^^o is highly enhanced in proximity of 
the transition II— >I. The result, confirmed by numerical analy- 
sis, thus indicates the presence of a divergent behavior at large 
D. This is a quantitative confirmation of the peculiar multipar- 
tite character of the transition II— T: at that transition, though 



We investigated in momentum representation the structure 
of correlations of the ground state of an extended ID Hubbard 
model, in relation to the undergoing QPTs. At variance with 
the standard real-space approach, the k-space representation 
allows in this case to explicitly derive the reduced density ma- 
trices of arbitrary subsystems, and thus to evaluate rigorously 
all bipartite and multipartite correlations present in the sys- 
tem. 

A first interesting result is directly related with the presence 
of single electrons with given spin and momentum; the conse- 
quence of the presence of these particles is the limitation of the 
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portion of momentum space available for the implementation 
of the eta-paring mechanism. We show that this fact has vari- 
ous important implications in the structure and the behaviour 
of the correlations present in the ground state. In particular, 
we show that the phase space of the model can be character- 
ized by the presence of iso-correlation curves, i.e. curves in 
the phase space along which the value and the structure of the 
correlations does not change. 

Due to the presence of the eta pairing mechanism, at a lo- 
cal level the existing network of correlations can be described 
in terms of two elementary building blocks: the generic sin- 
gle mode kj and the paired modes (kj, —kj). The identi- 
fication of such elementary subsystems and the comparison 
of their (quantum) correlations with the results obtained in 
the direct lattice picture allows one to unveil the origin of 
the off-diagonal long-range order (ODLRO), whose value can 
be functionally expressed in terms of quantum correlations 
between two paired modes, i.e., negativity . As a conse- 
quence, we also show that the mutual information between 
two sites can be directly linked with the negativity between 
paired modes. 

The identification of the above elementary subsystems has 
also allowed for the classification of the QPTs in terms of 
two-pair or multi-pair correlations, depending on whether par- 
tial derivatives of quantum mutual information and Von Neu- 
mann entropy show the same singularity or not at a given QPT 



These singularities can be related to those shown by second 
partial derivatives of the ground-state energy and, when possi- 
ble, to the the singularities of measures of correlations derived 
in the direct lattice picture. 

Finally, we have computed a proper direct measure of mul- 
tipartite entanglement, the Q-generalization of the Meyer- 
Wallach measure 20 - 22 The latter requires the evaluation of the 
reduced density matrices relative to subsystems constituted by 
blocks of modes of arbitrary size D. We have shown that, 
besides the different behaviour of the measure depending on 
the presence of paired/unpaired modes in the blocks 19 , as D 
grows the partial derivative of Q develops a singularity in cor- 
respondence with those QPTs at which the dominant corre- 
lations were conjectured to be of a multipartite nature. On 
the other hand in correspondence with the QPTs classified as 
"two-pair", by increasing D the partial derivative of Q be- 
comes smooth. 

The comparison between the structure of correlations in 
momentum space and their direct lattice counterpart turned 
out to be a useful method that allowed to unveil interesting re- 
lations between the correlations in the different pictures and to 
accurately characterize the ground state phase space and the 
quantum phase transitions occurring in a strongly correlated 
electron system. We expect that such method can be fruitfully 
extended to other many body complex systems. 
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